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ABSTRACT 

This paper presents a new sample of BL Lac objects selected from a deep (30 
mJy) radio survey of flat spectrum radio sources (the CLASS blazar survey, henceforth 
CBS). The sample is one of the largest well defined samples in the low power regime 
with a total of 130 sources of which 55 satisfy the 'classical' optical BL Lac selection 
criteria, and the rest have indistinguishable radio properties. The primary goal of this 
study is to establish the Radio Luminosity Function (RLF) on firm statistical ground 
at low radio luminosities where previous samples have not been able to investigate. 
The gain of taking a peek at lower powers is the possibility to search for the flattening 
of the LF which is a feature predicted by the beaming model but which has remained 
elusive to observational confirmation. In this study we extend for the first time the 
BL Lac RLF down to very low radio powers ~ 10 22 W/Hz, ie, two orders of magnitude 
below the RLF currently available in the literature. In the process we confirm the 
importance of adopting a broader, and more physically meaningful set of classification 
criteria to avoid the systematic missing of low luminosity BL Lacs. Thanks to the 
good statistics we confirm the existence of weak but significant positive cosmological 
evolution for the BL Lac population, and we detect, for the first time the flattening of 
the RLF at L ~ 10 25 W/Hz in agreement with the predictions of the beaming model. 
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1 INTRODUCTION 

BL Lacs are flat radio spectrum AGN which show high vari- 
ability in both flux and polarization and whose Spectral 
Energy Distributions (SED) are dominated by non-thermal 
processes. In these features BL Lacs are not alone. In fact, 
as early as the first conference dedicated to them as a class 
(Wolfe, 1978), the new term 'blazar' emerged, thus group- 
ing BL Lacs and Flat Spectrum Radio Quasar (FSRQs). 
Though this early amalgamation aimed at a unified model 
that would explain both types of sources, BL Lacs have nev- 
ertheless remained an elusive type of AGN throughout. For 
example, due to the lack of strong emission features in their 
optical spectrum, several BL Lacs are still without a red- 
shift measurement, despite continued efforts to obtain them 
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(see for instance, Meisner & Romani, 2010; Rau et al., 2012; 
Sbarufatti et al., 2009). Another aspect of the enigmatic na- 
ture of BL Lacs is the fact that there have been claims that 
a fraction of these sources suffers negative cosmological evo- 
lution (see for instance Morris et al., 1991; Wolter et al., 
1994; Bade et al, 1998; Giommi, Menna & Padovani, 1999; 
Rector et al., 2000; Giommi et al., 2012) something that is 
in marked contrast with the other blazars. 

Though there is good observational evidence to support 
the generally accepted view that the 'blazar phenomenon' is 
the consequence of observing radio galaxies down their rel- 
ativistic jets (see for instance Urry & Padovani 1995 and 
references therein; Capetti et al. 2002; Trussoni et al. 2003; 
Giommi et al., 2012, there are still detailed issues to be re- 
solved concerning this 'zeroth order unification'. In particu- 
lar it is important to keep in mind that though progress has 
been made in the analysis and high frequency observations 
of this type of AGN, it is still true that there is a lack of 
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well defined samples that allow the estimate of statistically 
important parameters such as the luminosity function (LF) 
and cosmological evolution, especially in the low luminosity 
regime. For instance, even though there are a few thousand 
of BL Lacs known (Massaro et al., 2011 mentions roughly 
1200 with the last version of the BZ catalogue), the num- 
ber of sources belonging to well defined samples from which 
statistical results can be obtained is much smaller than this. 
Furthermore, until recently these samples did not sample 
the population homogeneously in their radio luminosity. In 
fact, the high radio luminosity end of the population was al- 
most exclusively sampled via radio surveys, specifically the 
Uy sample (Stickel et al., 1991), whereas the lower radio 
luminosity samples were either X-ray selected or obtained 
through a combination of X-ray and radio data (see for 
instance Laurent-Muehleisen et al., 1998 and 1999 for the 
RGB; Perlman et al., 1998 for the DXRBS; Caccianiga et 
al., 2002b for REX). Since both types of surveys are sensitive 
to different parts of the SED, studying the entirety of blazar 
phenomenon could be compromised by selection biases due 
to the different ways the objects in low and high luminosity 
regime were selected (eg. see the discussion in Caccianiga & 
Marcha, 2004). 

Clearly the need to have larger samples that provide 
blazar data across the radio luminosity range is acknowl- 
edged (see for instance Giommi et al., 2012 for a discussion). 
The problem however resides in the very nature of BL Lacs 
since the lack of strong emission features in their optical 
spectra means that BL Lac samples are plagued with selec- 
tion effects. There are two main reasons for this: (1) the ini- 
tial BL Lac optical criteria required the EW of the strongest 
emission lines to be below 5A, something that is not only 
arbitrary, but also without much physical meaning, and (2) 
the fact that how weak a line is to be recognized before be- 
ing dismissed depends on a combination of instrumental and 
intrinsic source parameters. Several studies investigated the 
importance of the biases introduced in BL Lac samples due 
to these type of selection effects, and updated selection crite- 
ria considered more physically meaningful were introduced 
(Browne & Marcha, 1993, Marcha & Browne, 1995, Marcha 
et al., 1996, Landt et al., 2002, 2004). 

Making use of the expanded selection criteria new BL 
Lac samples have been produced (see for instance Giommi 
et al. 2012, for a discussion of different samples). The ad- 
vantage of these samples has been to broaden the parameter 
space within which the BL Lac population could be studied. 
Nevertheless, important parameters such as the Radio Lu- 
minosity Function (RLF) of BL Lacs has only been possible 
in two previous samples, both of which contain a small num- 
ber of sources - 34 sources in the case of the Uy sample of 
Stickel et al., 1991, and 24 objects in the case of the DXRBS 
from Padovani et al., 2007 - and reaching radio luminosities 
of around 10 25 W/Hz in the first instance, and 10 24 W/Hz in 
the second. 

In the current work we present a large and well-defined 
radio selected sample of low power BL Lacs. We use it to 
investigate the cosmological properties of this class of AGN, 
and to derive the Radio Luminosity Function (RLF) to radio 
luminosities as low as 10 22 WHz~ 1 . The paper is organised 
as follows: The next Section discusses the selection of the 
CLASS sample of blazars and how we used the available 
observational parameters to improve the BL Lac classifi- 



cation at low radio luminosities. Section 3 is dedicated to 
assessing cosmological evolution, whereas Section 4 sees the 
analysis of the RLF for the samples identified in the pre- 
vious sections. Finally, in Section 5, we discuss the results 
and conclusions are drawn. Throughout the paper we have 
assumed S„ <x v~ a " , with a a = 1.0 for the optical spec- 
tral index, and Ho = 71 Km s _1 Mpc _1 . However, when it 
came to comparing our results to those of previous studies 
we adopted Ho = 50 Km s _1 Mpc _1 , as some of the reference 
quantities have been determined with this value. We have 
also adopted an 'old fashioned' cosmology £l m — 0, fl\ = 
as these will have little effect on the low redshift regime we 
are primarily investigating. 



2 THE CLASS SAMPLE OF BLAZARS 

The CLASS blazar sample (CBS) consists of 302 flat spec- 
trum radio sources selected according to the following crite- 
ria: 

(i) 35° < S < 75° 

(ii) 1 6" | >20° 

(iii) S 5 > 30 mJy 

(iv) flat spectrum, i.e. a\% <0.5 (S„ tx v~ a ) 

(v) red magnitude < 17. E0 

The details of the selection, as well as the radio proper- 
ties of the sample are discussed in Marcha et al. (2001) and 
the reader is referred to that work for more detailed infor- 
mation. Further work on radio imaging and multi-frequency 
data of the sample enabled a refinement of the blazar classifi- 
cation to 244 (decreased from 302) of the objects in the CBS. 
This analysis is discussed in Caccianiga & Marcha (2004). 

Optical classification has been collected, either from 
the literature or from specific observations for about 91% 
of the objects in the sample. A discussion on the spectro- 
scopic classification is presented in Caccianiga et al. (2002a) 
but essentially, objects have been broadly divided into three 
categories: Type if their spectra show only weak or no 
emission lines (see next Section); type 1 if the spectra show 
broad emission features, and finally, type 2 if the spectra 
only show narrow emission lines. 

The focus of the current work is on BL Lacs, hence on 
sources whose spectrum is devoid of strong emission lines. 
However, at lower radio luminosities isolating these sources 
is not as trivial as it seems when dealing with more power- 
ful objects. In particular, two problems make this task more 
difficult: (1) the lack of an absolute measure of 'weak' and 
'strong', and (2) the difficulty of recognising weak BL Lac 
nuclei in the core of luminous host galaxies. In the absence of 
physically meaningful cut-off criteria different authors have 
made different attempts to minimise the effects of these two 
problems on the selection of a given sample. The approach 
followed in the present work is to investigate different sub- 
samples resulting from applying variations on the selection 
criteria that have been used in the case of previous samples 
of BL Lacs, and investigate how that affects the results and 
interpretation of their cosmological properties. A Table (??) 

1 Such a limit on the magnitude of the sources ensures that we 
target low redshift and low radio power sources 
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containing the relevant parameters of the sample used in the 
present work can be found in appendix A. 

2.1 The 'classical' BL Lacs 

Traditionally, a source is classified as a BL Lac if its optical 
spectrum is featureless or shows only weak emission lines. 
A working upper limit on the Equivalent Width (EW) of 
the strongest emission line has been set to 5A in order to 
separate BL Lacs from other radio-loud AGN. An additional 
criterion for the classification of BL Lacs in the EMSS survey 
was introduced by Stocke et al. (1991). It consisted in requir- 
ing that enough non-thermal continuum from the AGN was 
present in addition to the starlight of the host galaxy so that 
the observed contrast across the Call H&K discontinuity at 
4000A(A) was reduced to < 25%. However, the imposition 
of such a limit on the relative strength of the optical flux 
of the AGN to the starlight of the host galaxy is somewhat 
arbitrary and unlikely to be related to the intrinsic prop- 
erties of the active nucleus. In fact, the measured contrast 
is inversely proportional to the strength of the non-thermal 
contribution from the AGN, which means that as the AGN 
gets weaker, the measured contrast will get larger. Hence, 
the deeper the survey, the more important is the choice of 
contrast used in the selection of BL Lacs. A too limiting 
value for this quantity can introduce significant incomplete- 
ness in the final sample as is widely discussed in Browne & 
Marcha (1993) and Marcha & Browne (1995). 

It was first proposed by Marcha et al. (1996) while 
studying a sample of flat radio spectrum sources with flux 
limit of 200mJy that the limiting contrast for BL Lac defi- 
nition should be extended to 40% in order to avoid missing 
weak nuclei in the core of bright ellipticals. This was sup- 
ported by the lack of strong evidence to separate sources 
below and above the 25% threshold in a sample of optically 
bright, low luminosity flat radio spectrum objects. Since 
then, various works have lent further support for this ex- 
panded classification criterion, and a limiting contrast of 
40% has now been widely used in the definition of new BL 
Lac samples (Laurent-Muehleisen et al. 1998; Caccianiga et 
al. 1999, Perlman et al., 1996: Beckmann et al., 2003; Landt 
et al. 2002, Padovani et al., 2007). 

2.2 The Type sources 

Even though the inclusion of sources with A < 40% in the 
BL Lac classification prevents significant incompleteness of 
radio samples down to a flux limit of ~200mJy, it is not 
guaranteed that this limit will prevent weak BL Lac nuclei 
from being detected in much deeper surveys like the CBS 
where most of the selected low power BL Lacs are likely to be 
indistinguishable from normal "passive" elliptical galaxies. 
For this reason, we chose to study also all the sources in 
the CBS that are spectroscopically classified as Type 0, i.e, 
sources showing an optical spectrum with weak (EW<5A) 
or no emission lines, and no restriction on A 

There are 93 sources that fall into this category: the 
55 'classical' BL Lacs plus 38 sources named PEGs (Pas- 
sive Elliptical Galaxies) which apart from measured con- 
trasts above 40%, are otherwise indistinguishable from BL 
Lacs. There are good reasons for doing this as previous data 



such as the analysis of SEDs (Anton et al. 2004; Anton 
& Browne, 2005), radio imaging and polarization studies 
(Dennett-Thorpe & Marcha, 2000; Bondi et al., 2004), where 
at least in the case of some PEGs, it was not possible to dis- 
tinguish the morphology and polarization properties from 
those of 'classical' BL Lacs. This may indeed be the best 
observational evidence yet to support the existence of a BL 
Lac nucleus in the core of at least some PEGs. In the case of 
the CBS it is more difficult to obtains such a direct confir- 
mation of the presence of a BL Lac nucleus in PEGs due to 
the weakness of the selected sources. However, we have pre- 
sented evidence, based on high resolution VLA maps, that 
this is likely to be true also in this sample (see Caccianiga 
& Marcha, 2004). 

2.3 The Weak Emission Line sample 

The optical criteria separating BL Lac objects from quasars 
is known to be somewhat arbitrary, and in the two previous 
sections we discussed ways to prevent incompleteness in BL 
Lac samples due to this arbitrariness especially when one 
has the advantage of starting from a radio survey. Other at- 
tempts to address the optical identification of BL Lacs have 
however been published in the literature, and it is important 
to investigate how they may compare to those discussed in 
the two previous Sections. 

In 2004, Landt et al. adopted a diagnostic plot based on 
the equivalent width of two narrow emission lines, namely 
the [OII]A3727A and the [OIIIJA5007A. According to these 
authors, this diagnostic plot is able to separate AGN that 
are intrinsically "weak lined" (WLAGN) from the "strong 
lined" ones (SLAGN) independently to their orientation. 

For consistency with what is presented in the literature 
we want to apply also this alternative classification scheme 
to the CBS sample to produce the RLF of the population 
of WLAGN which should be considered an extension of the 
population of type based on a more "physical" classifica- 
tion. 

Only for ~ half of the 244 CBS sources classified as 
"blazar", based on radio data (see Caccianiga & Marcha 
2004 for details), do we have a spectrum covering both 
the [OII]A3727A and [OIII]A5007A emission lines used in 
the work of Landt & co-workers. The sources for which we 
have not covered one of the two lines (or both) are typ- 
ically high z objects (z>0.8). However, the choice of the 
classification scheme is really critical for blazars with a rel- 
atively low luminosity (and, hence, low redshift), where 
the possible effects due to the orientation may be rele- 
vant. For powerful blazars the Landt et al. (2004) converges 
to the usual classification based on the EW=5A thresh- 
old. An indicator of the AGN optical luminosity is the 
value of the 4000A break (A): sources with A <20% are 
usually powerful blazars (typically with Psghz >10 25 W 
Hz -1 ) for which the Landt et al. (2004) classification in 
WLAGN /SLAGN is almost equivalent to the Type 0/Type 1 
classification that we have adopted in the CBS sample. For 
these sources with A <20%, we have thus simply "con- 
verted" our classification into the WLAGN/SLAGN classi- 
fication by assuming type 1=SLAGN and type 0= WLAGN. 
For the sources with A >20%, instead, we have applied 
the [OII]A3727A/[OIII]A5007A classification plot proposed 
by Landt et al. (2004) to distinguish between WLAGN and 
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Figure 1. The classification plot proposed by Landt et al. (2004) 
applied to the CBS blazars with large values of 4000A break 
(A >20%). The solid line is the dividing line proposed by Landt 
ct al. (2004) to separate between WLAGN and SLAGN. Symbols 
indicate the former classification discussed in Caccianiga et al. 
(2002a): filled circles = type 0, open circles = type 1, stars = 
type 2 



SLAGN. Since these sources are typically at low redshifts 
(z<0.3) both emission lines are covered by our spectra in 
the large majority of cases (~80%). For the remaining ob- 
jects the missing emission line is typically the [OII]A3727A 
because the available spectrum does not cover the blue part. 
As discussed in Land et al. (2004) the [OIIL1A5007A alone 
can be used in these cases to provide a tentative classifica- 
tion: if it is large (>20A), then the source is likely a SLAGN 
while, if it is small (<3A), the source is likely a WLAGN. In 
the intermediate cases the ambiguity can be removed only 
by measuring the [OII]A3727A equivalent width, so we con- 
sider these sources as "unclassified". Only 17 blazars, out 
of 244 (7%), belong to this group of unclassified sources. In 
Fig. Q] we report the classification plot for the CBS sources 
with A >20% and for which both emission lines are covered. 
Most of these objects fall in the WLAGN region. 

In total, the sample of WLAGN includes 129 objects. 
This sample contains all but one of the type objects, de- 
scribed in the previous section, plus some low-power sources 
previously classified as type 1 or type 2 objects. We have 
also decided to include the two Type objects whose upper 
limits on the EW of the two lines he above the dividing line 
in Figure □} as WLAGN. 



3 ASSESSING COSMOLOGICAL EVOLUTION 

The evolutionary properties in this work are investigated via 
the 'V/Vmax test'. First introduced by M. Schmidt (1968) 
to study the space distribution of quasars the test was later 



generalised by Avni & Bachall (1980) for the case where the 
sample selection involved more than one survey. Basically, 
the test consists in assessing the ratio between the volume 
V , corresponding to the volume "enclosed" by the redshift of 
the source, and the maximum volume (Vmax) within which 
the source can still be found satisfying the selection criteria 
of the survey. In the specific case of the CBS there are two 
flux limits that need to be considered, one in the radio and 
one in the optical, and consequently, the Vmax to be used 
will be the one corresponding to the most limiting of the 
two. 

For non-evolving populations, the ratios V/V m ax com- 
puted for each source should be uniformly distributed be- 
tween and 1 with a mean value < V/Vmax >= 0.5 and sta- 
tistical errors being given by o = l/yf(12N) (where N is the 
number of sources in the sample). If, on the other hand, the 
distribution of sources is non-uniform, < V/Vmax >^ 0.5, 
and the population is said to be evolving: positively if the 
mean is larger than 0.5, or negatively if it is smaller than 
0.5. It is however, important to emphasise the fact that a 
deviation from a uniform distribution of V/Vmax can also 
originate, not due to evolution but rather from unrecognised 
selection effects in the sample. To distinguish between the 
two situations is particularly relevant in the case of BL Lacs, 
where primarily due to the properties of their optical spec- 
tra, weak nuclei can easily be missed in flux limited samples. 

Before proceeding with the analysis, there is one issue 
that must be addressed: how to deal with the z-less BL 
Lacs. There are 16 such sources which is a significant fraction 
of the 55 BL Lac sample. We have explored a number of 
possibilities for attributing a redshift to each of the 16 z-less 
BL Lacs. We started out by considering just the subgroup 
of 39 BL Lacs with z measurements and then the following 
possibilities: (1) all z-less sources were assigned a z = 0.3 
(the mean of the BL Lacs with measured redshift, (ii) a 
z = 0.8 for all the z-less BL Lacs, and finally (iii) a random z 
was found for the 16 z-less sourcesQ The results can be found 
on Table [1] Essentially, there is minimal impact on the < 
V/Vmax > as a result of assigning different values of z to the 
z-less sources. Consequently, for the remainder of this work 
we have considered the option where the 16 z-less BL Lacs 
have been assigned a random value for their redshift (drawn 
from the LI interval as defined in the footnote). We believe 
that using a random distribution for the z-less sources is, not 
only a more realistic approach than using the mean value of 
the sample for all the objects without redshift and one that 
is likely to have the least fictitious effects on the estimated 



2 This was done by allowing the sources to have radio luminosi- 
ties in the range typical of that observed for other radio selected 
BL Lacs, e.g. LI: 10 24 < L r < 3 X 10 28 W/Hz with an extra con- 
dition of z > 0.1 since below this value (assuming a typical radio- 
to-optical ratio for BL Lacs), we should detect the host galaxy 
in the optical image. We have taken the exercise a step further 
and actually investigated two further luminosity intervals from 
where the random z distribution of the z-less objects were drawn: 
L0:10 23 < L r < 3xl0 28 , and L2 :10 25 < L r < 3xl0 28 W/Hz. We 
found that apart from the highest z achieved (z m ax ~ 3 for L0, 
to Zmax ~ 1.5 for L2), it had little influence on the < V/Vmax > 
test, e.g. mean and distribution. We have used the interval LI 
as a default but will discuss the implications of using the other 
luminosity intervals when analysing the LF. 
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RLF, but also one that is in agreement with a recent work 
by Rau et al. (2012) which suggest that some z-less BL Lacs 
have photometric redshifts as high as z ~ 2.. 

We have analysed 4 subsamples separately, 3 from the 
CBS - BL Lacs, TypeO, WLAGN - plus the Uy BL Lac 
sample of Stickel et al., 1991 as it provides a validity test on 
the analysis. 

3.1 Do BL Lacs evolve? 

The cosmological evolution of BL Lacs has been a matter 
of debate over the last two decades and it is still considered 
as an open problem. At the origin of the debate are the op- 
posite results obtained when using the only 2 statistically 
complete samples of BL Lacs available in the early 90s, one 
selected at radio frequencies (1 Jy sample; Stickel et al. 1991) 
and the other in the X-ray (the EMSS-C sample; Morris et 
al. 1991). In the case of the radio selected sample the result 
was consistent with a population that suffered no evolution, 
or a mildly positive one, whereas in the case of the X-ray 
selected sample the result was clearly indicative of a popu- 
lation that evolved negatively. The evolutionary picture of 
BL Lacs became even less clear when later results based on 
ROSAT data (Bade et al. 1998; Giommi, Menna & Padovani 
1999; Beckmann et al., 2003; Giommi et al. 2001) seemed to 
confirm the presence of negative evolution for a subset of the 
BL Lac population, i.e. for the High energy peaked BL Lacs 
(HBL) , while other sets of data, namely the DXRBS and the 
REX surveys found no evidence for evolution (Caccianiga et 
al., 2002b; Padovani et al., 2007). 

More recently, Giommi et al., (2012), suggest that BL 
Lacs are essentially of two types: those that are beamed 
cores of FRII radio galaxies and which should be grouped 
with FSRQs, thus sharing their strong cosmological evolu- 
tion, and those that show no evolution because they are 
beamed versions of FRIs. The authors performed extensive 
numerical simulations and found that if they consider only 
the most extreme HBLs, then an apparent negative cosmo- 
logical evolution appears as a result of a combination of dif- 
ferent selection effects. Can the CBS sample help shed any 
light onto this matter? 

Table [5] shows the values obtained for the < V/Vmax > 
test for each of the samples mentioned before. It also shows 
the results of modelling the evolution with a simple expo- 
nential form exp(C x r(z)) where t(z) is the look-back time 
in units of H . The values of the K-S test for the distribu- 
tion of V/Vmax before and after the evolution model are also 
shown. 

The results can be summarised as follows: 

(i) We used the Uy sample of BL Lacs of Stickel et al. 
(1991) for a reliability check on our analysis and we find sim- 
ilar values for all the parameters (the same < V/Vmax >= 
0.6 ± 0.05, and C = 3 vs. C = 3.1 cited in Stickel & co- 
workers). More interestingly, the sample shows the weakest 
statistical need to invoke evolution. With a 2a detection and 
a V/Vmax distribution that cannot be dismissed as incom- 
patible with a uniform one, it is the sample where a LE 
model was fitted with the lowest exponential C parameter. 

(ii) The CBS samples of 'classical' BL Lacs, TypeO and 
WLAGN show similar < V/Vmax > values and V/Vmax 
distributions although BL Lacs showed the highest, and 



WLAGN the lowest deviation from the mean value of 0.5. In 
both the CBS 'classical' BL Lac and TypeO samples the de- 
viation is now at 3a whereas for the WLAGN the deviation 
is only detected at 2a. Also, the V/Vmax distributions are 
now inconsistent with being uniform, contrary to the case 
in the Uy sample. 

(iii) Modelling the evolution yielded similar values for 
an LE exponential model although the BL Lacs required 
the highest evolution parameter and the WLAGN the low- 
est. We note however, that the evolution parameter C that 
achieves < V/Vmax >= 0.5 within la is pretty similar in the 
3 cases; 

(iv) How do these values compare with other works? The 
values of < V/Vmax > obtained for the CBS samples of BL 
Lacs, WLAGN and TypeO are comparable to the one found 
for the Uy sample of BL Lacs but whereas in the latter the 
deviation was only detected at 2a, in the case of the CBS 
samples the numbers allow a 2 or 3a confidence (depend- 
ing on the sample in question) in the claim of cosmological 
evolution. Furthermore, the V/V max distributions are not 
consistent with a uniform one, in contrast to what happens 
in the Uy sample. We also found that the evolution param- 
eter values required to bring the < V/Vmax > down to 0.5 
were larger in the CBS samples (C ~ 5 — 6) when compared 
to that obtained for the Uy sample (C ~ 3). 

Apart from the Uy sample of Stickel et al. (1991), the 
closest comparative study is the one carried out by Padovani 
et al., 2007 on the DXRBS blazar sample. In this study 
Padovani and co-workers find that the 24 BL Lacs that con- 
stitute the sample show < V/Vmax >= 0.57 ± 0.06, a value 
consistent with the one found for the CBS BL Lacs. We 
therefore conclude that the BL Lac population suffers weak 
(when compared to FSRQs) but significant positive evolu- 
tion, a result that was difficult to obtain due to the small 
statistics of previous samples. 



4 THE RADIO LUMINOSITY FUNCTION 

The differential luminosity function is determined according 
to the following equation: 

N 

where V m ax is the maximum volume, A corresponds to 
surveyed area, AL the width of the luminosity bin, and N 
is the number of objects contained in the bin (L,L+AL). 
As discussed previously in Section 3, Vmax is the smallest 
maximum volume between those computed using the radio 
and the optical limits of the sample. In this way the presence 
of the optical magnitude cutoff is automatically taken into 
account also in the derivation of the RLF. 

4.1 The analysis of the RLF 

The analysis of the LF holds vital statistical information 
about the population being studied. For example, some au- 
thors claim that the slope of the RLF of radio galaxies can 
give an idea of the rate at which AGN producing jets ap- 
pear in the universe (Kaiser & Best, 2007). In the particular 
case of BL Lacs, the underlying assumption is that these 
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BL Lacs w/ z z-less= 0.3 z-less= 0.8 z-less=random 
< V/Vmax > 0.59 0.62 0.63 0.62 

ly/(12N) ±0.05 ±0.04 ±0.04 ±0.04 

Table 1. Mean V/V m ax values for the BL Lac CBS sample for different redshift options for the z-less sources. 



Sample 


< V/Vmax > 


K-S (%) 


evol. model 


C 


< V/Vmax > 


K-S (%) 


BL Lacs 


0.62 (±0.04) 


«1% 


LE(exp) 


6.1i: 5 


0.50 


~60% 


Type 


0.59 (±0.03) 


«1% 


LE(exp) 




0.50 


70% 


WLAGN 


0.57 (±0.03) 


<1% 


LE(exp) 


5.±i: 


0.50 


99% 



Table 2. Summary of the parameters for the evolutionary model in the case of the different samples. Columns are: (1) Sample; (2) 
< V/Vmax > obtained and its respective error; (3) K-S test probability for the < V/V m ax > (see text for details); (4) type of evolutionary 
model used, (5) Evolutionary parameter where the upper and lower limits indicate the C which still falls within < V/Vmax >= 0.5 ± la; 
(6) K-S test probability for the < V/Vmax > distribution after the evolution was taken into account. 



sources consist of beamed up cores of low luminosity radio 
galaxies. In such a scenario the LF of the two populations 
should be connected. Basically, if the LF of the parent pop- 
ulation (unbeamed) consists of a single power-law, the effect 
of beaming produces a break in the LF at a luminosity that 
depends on the Doppler factor S — [F(l — /3cos#)] _1 , where 
r = (1 — /3 2 ) -1 / 2 is the Lorentz factor, /3 the jet velocity 
in units of the velocity of light, and O the angle between 
the jet and the line of sight (l.o.s). Detailed analysis of the 
effect of beaming on the LF has been carried out by several 
authors (see for example Urry & Shafer, 1984 and Urry & 
Padovani 1991 and 1995 for thorough discussions). 

In the present work the RLF for the CBS will be stud- 
ied under two scenarios: evolution (with an exponential LE 
model), and no-evolution. For each case, the LF was fitted 
by: 



$(L) = KL a 

for a single power-law, or 
$ 



*(L) 



(2) 



(3) 



for a broken slope at a given luminosity Lbreak- The 
results are shown in Table and Figures [2] [3] and [4] 
The results can be summarised as follows: 

(i) For the 'classical' BL Lacs of the CBS we find that 
the reasonable fits to the shape of the RLF are those of a 
single power law (see Figure [4] and Table The slope in 
the evolution scenario is consistent with the steepest value 
found for the Uy sample of BL Lacs (Stickel et al., 1991 
found slopes between 2.4 and 2.7 depending on the binning), 
and those found for the DXRBS (slope of ~2.3 for a bin 
size of AlogP = 0.6). The evolution parameters are quite 
different in the three samples (C = 6 in the case of the CBS 
and C = 3 and C ~ 2.2 for the 1 Jy and DXRBS samples, 
respectively). It is important to note that just as in the case 
of the DXRBS, the CBS 'classical' BL Lac sample achieves 
radio luminosities at least a factor of 10 below those of the 
Uy BL Lac sample. The total number density found for the 
CBS 'classical' BL Lac sample from the integral LF for the 
evolution scenario is between those values found for the Uy 



and the DXBRS. Specifically, for the same cosmology (fl m = 
0, fl x = 0, and H = 50 Km s^Mpc" 1 ) Stickel et al. (1991) 
found a value of N ~ 40 Gpc -3 for the interval of luminosity 
interval 6 x 10 24 W/Hz - 3 x 10 27 W/Hz, while Padovani & 
co-workers found the total number density of the DXRBS 
BL Lacs to be N ~ 310 Gpc -3 for the the luminosity interval 
of 10 24 W/Hz to 6xl0 2B W/Hz. For the CBS 'classical' BL 
Lac sample this number was found to be N - 90 Gpc -3 for 
luminosity range 10 24 W/Hz to 6xl0 26 W/Hz. 

(ii) For Type and WLAGN we find that the RLF is best 
fit by a broken power-law (see Figure [4] and Table [3jl - It is 
important to emphasize that the RLF for the TypeO and 
the WLAGN show good overall agreement in the region of 
overlap with the Uy sample (i.e., luminosities roughly be- 
tween 10 25 and 3 x 10 27 W/Hz - see Figure but that their 
luminosity range extend to a factor of ~ 1000 times below. 
We note that such an agreement warrants credibility to the 
methodology used to derive the statistical properties of the 
sample. This is particularly relevant for the high luminos- 
ity end of the LF which is the region where the CBS starts 
loosing coverage, and one may worry about the sample not 
being representative, either because of the magnitude limit 
in the selection criteria, or because of the assumptions made 
for the sources without redshift measurements. The fact that 
the CBS samples reach all but the very last luminosity bin 
of the 1 Jy BL Lac luminosity function without any real 
discrepancy between the two, gives good guarantee that any 
missed sources at the high end of the luminosity range would 
have little impact on the statistical analysis of the samples. 
The total number density (for the evolution scenario and 
Ho = 50 Km s~ Mpc - cosmology considered before) de- 
rived from the integral LF for the luminosity range 3 x 10 22 
- 3 x 10 26 W/Hz was found to be N ~ 1.4 x 10 4 Gpc" 3 and 
N ~ 5 x 10 4 Gp" 3 for the TypeO and WLAGN, respectively. 

(iii) Taking the TypeO /WLAGN we find for the first time 
evidence of the break in the LF according to the predic- 
tions of the Unified Schemes. Indeed, the broken power law 
fit is statistically preferred in respect to the simple power 
law fit as demonstrated by the F-test (probabilities below 
5-10%). More specifically, starting from the RLF of FRI 
galaxies, Urry, Padovani & Stickel (1991) predict that the 
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Figure 2. Estimated RLF in the evolution scenario. From Top to bottom: BL Lacs, TypeO, WLAGN. The best fits are shown by the 
continuous lines. (Hq = 71 Km s — 1 Mpc — 1 ) 
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Figure 3. RLF for the TypeO and WLAGN samples (filled circles) against the 1 Jy sample (open squares) of BL Lacs in the no evolution 
scenario. (Hq = 50 Km s~ 1 Mpc _1 ) 
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Sample fit (no cvol.) x 2 /Dof (evol) x 2 /Dof 



BL Lacs 
p-law 

broken p-law 



a 
b 

Log L b (W/Hz) 
Log Lb range 



-2.09 



+0.08 
-0.08 



4/7 



C=6 



-2.31 



+0.12 



7/4 



Type 
p-law 

broken p-law 



-2.02 1 



2.49" 1 



1.66" 1 



Log L b (W/Hz) 25.09±™4 
Log L b range [24.75,25.21] 



23/12 
10/10 



C=5.5 

+0.07 



-2.15 



3.3V 



+0.12 



1.81 



25.0ll°jS 



[25.01, 25.67] 



37/9 
21/7 



WLAGN 
p-law 

broken p-law 



Log L 6 (W/Hz) 
Log Lj, range 



-2 07°' 04 
z ' u '-0.04 



,+0.47 



2.71] 
1.87j 
25.i 
[25.69,26.54] 



-0.27 

7 +0.10 

-0.12 

3 +0.63 
="-0.76 



25/10 
13/8 



C=5 



-2.29 



+0.05 
-0.06 



2.93 
1.93 
24.64 



+0.33 



0.24 

+0.15 

0.19 

-0.51 



-0.56 

[24.64, 24.99] 



24/8 
8/6 



Table 3. Table containing the parameters for the fitting of the LF. Columns are as follows: (1 and 2) - Sample and parameters for the 
type of fit, either power-law or broken power-law; (3 and 5) - slope of the LF, luminosity break, and the range of the luminosity break 
obtained by varying the luminosity interval from which the random z distribution was obtained, e.g L0-L2 as defined in Section 3.; (4 
and 6) - x 2 /degrees of freedom. 



beamed RLF has negative slopes of around 1.6 for lumi- 
nosities between 10 21 — 10 25 W/Hz, and 3.3 for the range 
4.7 x 10 27 - 10 29 W/Hz. It is remarkable that without any 
parameter fixing, these values are in such relatively good 
agreement with those found for our fits of the RLF for 
the TypeO /WLAGN (see Tableland Figure gj. Based on 
their simulations, Urry, Padovani & Stickel (1991) find a 
space density for BL Lacs (Ho = 50 Km s _1 Mpc _1 ) is 
N ~ 5.1 X 10 3 Gpc" 3 for sources with L > 5 x 10 20 W/Hz. 

(iv) The comparison between the RLF derived from the 
CBS samples and the prediction from the beaming model 
are also shown in Figure [4] Since the high luminosity end 
of the RLF is critically dependent on the adopted value for 
the evolutionary parameter, we have used the same param- 
eter C — 3 used by Urry Padovani (1995) to derive the 
beaming parameters. We have also applied the correction 
factor of 3/2 in the normalization discussed in Padovani et 
al. (2007) to account for the different classification criteria 
adopted. What we observe is that the RLF of the TypeO and 
WLAGN follow nicely the predicted (but never observed be- 
fore) shape of the LF, at least down to L ~ 10 23 W/Hz. A 
significant departure from the prediction is only observed in 
the low end of the luminosity range (L ~ 10 22 - 10 23 W/Hz). 
Such a feature can imply that either the parent population 



of BL Lacs was not sampled correctly at the low end of the 
luminosity range at the time of Urry & Padovani's work, or 
that the low luminosity end of the CBS RLF is contami- 
nated by "interlopers" (a situation already discussed in the 
Caccianiga & Marcha, 2004). Deeper observations of these 
very low-power sources (only about 5 objects in total) are 
thus required. In any case, the presence of the luminosity 
break (at much higher luminosities) is not affected by this 
low L tail. 

It is interesting to note that although BL Lac objects with 
radio luminosities in the range between 10 23 - 10 24 W/Hz arc 
predicted by the beaming model, no 'classical' BL Lacs are 
actually found, only TypeO or WLAGN. This, we believe 
is a convincing indication that the usual (classical) BL Lac 
definition, simply does not work in this range of lower radio 
luminosities, and that PEGs or WLAGN should be consid- 
ered as BL Lacs. 

(v) We have investigated how the assumption we made 
for the z-less sources affects the predicted RLF by fitting the 
LF each of the three random z distributions resulting from 
the three different L intervals mentioned in Section 3. The 
results showed that the effect of changing the distribution 
of the z-less sources was to shift the value of Lt rea k- For 
example, for the case of the Type sources we found that 
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Figure 4. Points indicate the observed LF de-evolved according to an exponential LE model with evolution parameter C = 3 against 
the predicted LF according to Urry, Padovani & Stickel (1991) for the TypeO, WLAGN and 'classical' BL Lac samples discussed in this 
work. {H = 50 Km s^Mpc" 1 ) 



the LogLbreak varied between 25.01 and 25.67 for the case of 
evolution and Ho = 50 Km s _1 Mpc _1 . For the same scenario 
but for the WLAGN sample we found that the variation of 
LogLbreak due to the range in the z-less distribution (e.g the 
luminosity intervals used to derive the redshift distribution 
for the z-less sources L0-L2 as defined in Section 3.) was 
between 24.64 and 24.99. 



5 SUMMARY AND CONCLUSIONS 

We have used a well defined sample of low power flat radio 
spectrum sources from the CLASS catalogue with the prime 
objective of studying the radio luminosity function of BL 
Lacs. This objective lead us to investigate not only the 'clas- 
sical' BL Lacs found in the CBS, but also their closest rela- 
tives, ie, the other flat radio spectrum sources where optical 
line emission is clearly weaker than that found in quasars. 
The present investigation builds on knowledge gained from 
the extensive work carried out on a pilot sample, the so 



called '200mJy sample' (Marcha et al., 1996) - see for in- 
stance Bondi et al., 2004, Dennett-Thorpe & Marcha, 2000, 
Anton et al., 2004 - and which supports the idea that at 
least some of the weak lined galaxies with radio properties 
similar to the classical BL Lacs, are indeed weak BL Lac 
nuclei which reside in the nucleus of bright ellipticals. 

Hence, the aim of the current study was twofold: to 
scrutinise the RLF of BL Lacs at low radio powers whilst 
investigating the sources that lie near the cut-off limits of 
'classical' BL Lac definition. Allowing for more physically 
meaningful selection criteria, and in accordance with previ- 
ous works, we analysed three subsamples from the CBS: the 
'classical' BL Lacs, the TypeO and the WLAGN. 

The conclusions can be summarised as below: 

(i) We present one of the deepest and largest samples of 
BL Lacs selected at radio frequencies. This gives the un- 
precedented opportunity to explore the BL Lac population 
to radio luminosities over 100 times lower than existing sam- 
ples. Furthermore, the work also greatly improves the statis- 
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tics for the samples of these types of sources (nearly a factor 
of 4). 

(ii) Due to the improved statistics we are able to detect 
a weak positive evolution (only marginally detected in pre- 
vious studies) with a 2-3a confidence. 

(iii) Reaching lower radio luminosities has enabled the 
detection of the change in the slope of the RLF predicted 
by the unified schemes at roughly 10 25 W/Hz, but never wit- 
nessed before. The agreement between the observed RLF for 
the TypeO and WLAGN sources of the CBS and the model 
predictions are remarkable considering that there was no 
fitting involved. 
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Table Al. Table containing the sources used in this work. Columns are as follows: (1) Name; (2) Position (J2000); (3) redshift; (4) 
Radio flux at 5 GHz (mjy); (5) R magnitude; (6) radio spectral index; (7) Classification: bl for 'classical' BL Lac, tO for Type and w 
for WLAGN. 
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The CLASS BL Lac sample: The Radio Luminosity Function 



Table Al. continue 



name 






Position 
(J2000) 
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(mjy) 
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